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1. INTRODUCTION 

Population growth [1] and migration of rural population towards cities have directly influenced the 
growth of vehicle population [2] in urban areas leading to congestion and increased travel time for commuters. 
Using public transit services such as buses can reduce the aforesaid problems of commuting in tier-2 cities (in 
countries like India) with limited mass transit options. To pull commuters toward public transit, it is necessary 
to provide reliable transit service [3]. Advanced passenger information system (APIS), a part of the intelligent 
transport system (ITS) [4] aims to provide dynamic information to passengers about the schedule, arrival, and 
travel times making the public transit service more reliable. Providing real-time information to passengers 
requires the integration of several static data such as spatial characteristics of the road network, schedule of 
buses, headway, and real-time data sources such as live location information of the buses, traffic flow, 
incidents, and congestion-related information. In cities with limited ITS infrastructure the aforementioned real- 
time data sources are rarely available but, recently the public transit buses in many cities are equipped with 
global positioning system (GPS) [5] that provides the dynamic location of the buses, speed, and heading. 

The travel time naturally varies depending on the space and time such as the day of the week, 
(weekday or weekend), time of the day, the land use pattern (LUP) of the route, the number of signalized and 
non-signalized intersections, and so on. The LUP is correlated to traffic density, commercial/business areas 
exhibit high traffic volumes leading to congestion and more travel time. Also, the stochastic behavior of the 
traffic affects the travel time of buses. Delays at intersections [6] due to queuing, non-optimized signals [7], 
and mixed traffic lanes lead to an excess travel time for buses. Additionally, the bus driver’s behavior [8] in 
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scenarios such as early start, on-time start, and delayed start of the trips will make the prediction process more 
complicated. Public transit travel time is predicted in two scales, long [9] and short-term [10]. If the travel time 
is forecasted more than 60 minutes ahead of the current time, it is a long-term, else a short-term prediction. 
Long-term prediction aid the decisions for operations planning of the buses [11], while short-term assists the 
passenger information systems [12], bus routing, fine-tuning schedules, and identifying bus bunching [13]. 

In the existing works, the bus travel times are predicted using manually collected data [14], GPS logs 
of the buses [15], automatic passenger counters [16], and mobile phone footprints [17]. GPS-based location 
data is the most common data source used for bus travel time prediction. From naive prediction models such 
as historic averages [18], statistical models like Kalman Filtering [18], time series models such as 
autoregressive integrated moving average (ARIMA), and seasonal autoregressive integrated moving average 
(SARIMA) [19]. Machine learning (ML) models such as artificial neural networks (ANN) to deep networks 
models [20], several authors have researched on the application of these models for short and long term travel 
time predictions. Considering high variances in the travel time authors in [21] have implemented the support 
vector regression (SVR) model in Chennai, India. Spatial and temporal SVR model is implemented to predict 
dynamic travel time using only GPS data of buses. Authors in [22] have used probe vehicle data, segment data, 
and weather data in Charlotte, North Carolina for short-term travel time prediction. Decision trees (DT), 
random forest regression (RFR), extreme gradient boosting (XGBoost), and long short-term memory (LSTM) 
models are compared and it is found that RFR outperforms other models. In [23] non-linear and linear models 
are compared for predicting the bus travel time using location data of the public transit buses and authors found 
that non-linear models DT, RFR, gradient boosting regression trees (GBRT) and k-nearest neighbors (KNN) 
outperformed the linear models; Linear regression, SVR, least absolute shrinkage and selection operator and 
ridge regression models, highlighting the non-linear behavior of travel times. The study was conducted in a 
tier-2 city (Tumakuru), India. A study [24] in Mumbai, India estimated the time of arrival of buses at bus stops 
and intersections using parametric hazard models based on Cox regression and accelerated failure time (AFT). 
Log-logistic and Weibull distributions are fitted to estimate the arrival time of buses at bus stops and 
intersections, and the authors found that the AFT model performed better with a 10% reduction in prediction 
variation. 

Though the prediction of travel time looks like a simple regression, the randomness of traffic density, 
weather conditions, congestion, and incidents demand sophisticated prediction models. Several researchers 
have developed short-term travel time prediction models using ML but few authors stress the necessity of 
hybrid models. However, most of the existing works are conducted in tier-1 cities with matured traffic 
infrastructure. But in tier-2 cities with limited infrastructure and data sources scenario, a dynamic model to 
predict travel time is yet to be explored. In this work, short-term travel times of buses are predicted using 
limited data sources, i.e., historical location data of buses, schedule, bus-stop information, and road geometry 
data. The objectives of the study are to predict/estimate the travel time of public transit buses at two spatial 
aggregation levels; route and segment levels in a tier-2 city with limited data sources as follows: 

a. Apply the GBRT machine learning model based on the historic data to predict the bus travel time. 

b. Propose the latest travel time estimation model (LTTEM) to estimate the dynamic bus travel time utilizing 
the travel time information of the preceding bus and historical data. 

c. Propose a hybrid travel time estimation model (HTTEM) that combines the results of the GBRT and 
LTTE models to improve the performance as compared to the individual models. 


2. DATA 

The location data of public transit buses of Tumakuru city during March 2021 is used for the modeling. 
The data includes the bus number, device-ID, date and time, latitude, longitude, speed, odometer, and location 
information. Route number 201: Tumakuru bus stand (TBS) to Kyathasandra (KYA) is chosen for model 
development and validation. The route under study is split into four segments based on the LUP. Segment | is 
from the TBS to Bhadramma choultry (BC) is a part of the central business district (CBD). Segment 2 is from 
BC to Shivakumara Swamiji circle (SSC) is an inner-city (IC) area, segment 3 from SSC to Batawadi (BW) is 
the inner suburban (ISU) region, and finally from BW to KYA segment 4, which is an outer suburban area 
(OSU). The map of the bus route is given in Figure 1, highlighting the bus stops and location of the split of 
segments.The route information is given in Table 1. Sample location data are presented in Table 2. The raw 
location data needs to be cleaned and pre-processed for analysis. For each trip start time, end time, running 
speed, length traversed, and total travel time is extracted from the data at both the route level and the segment 
level. 500 upstream trips on weekdays between mornings 7:00 to evening 20:00 are used in the study. The 
LUP, the number of intersections, and the segment number are augmented to trip aggregates. The weather data 
is not considered as the study location has stable weather, and the climatic changes rarely impact the travel 
times during the study period. 
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Figure 1. Map of route number 201 


Table 1. Route information 


Parameters Sl $2 $3 S4 Route 
Origin — destination TBS-BC  BC-SSC SC-BW BW-KYA _ TBS-KYA 
Length in kilometers 1.76 1 2.09 2.05 6.9 

Bus stops 3 2 3 1 9 
Signalized intersections 3 1 1 1 6 
Land use pattern CBD IC ISU OSU 


Table 2. Sample location data of route 201 


Vehicleregno Date_time Location Speed Odometer Latitude _ Longitude 

KA-06-F-0857 01-03-2021 07:29:58 TBS 13 12673.91 13.34319 = 77.09875 

KA-06-F-0857 01-03-2021 07:30:08 TBS 22 12673.96  13.34275  77.09859 
Near by Karnataka state road transport 

KA-06-F-0857 01-03-2021 07:30:18 corporation bus depot 23 12674.03 13.34216 =—77.09837 
Near by Karnataka state road transport 

KA-06-F-0857 01-03-2021 07:30:28 corporation bus depot 8 12674.07 13.34181 = 77.09828 

KA-06-F-0857 _ 01-03-2021 07:30:38 Near by 4148, PH colony 21 12674.11 _13.34145 7.09815 


3. METHODS 

Estimating the dynamic travel time needs the historical data of the trips, the up-to-the-minute location 
data of the buses, traffic on links, delays at intersections, and any incidents in the link. In most tier-2 cities, 
these data are rarely available, but developing models to estimate dynamic travel time is crucial in such 
locations. In this context, the travel time of buses is estimated using the limited data available by applying the 
methods discussed in the following subsections. 


3.1. Gradient boosting regression trees 

In a study conducted [23] it is identified that non-linear models are suitable for predicting travel time 
in the current study location. Therefore, travel time data are modeled using a non-linear ML model, GBRT. As 
the name suggests it is an ensemble greedy model that works based on boosting approach. GBRT works on the 
notion that the possible next model, when combined with the preceding weak model (regression tree) will 
reduce the error in overall prediction, and every proceeding model tries to reduce the errors of the collective 
boosted ensemble of all preceding models. It is a widely used model in predicting travel times. 


3.2. Latest travel time estimation model 

The travel time of buses is highly influenced by the real-time traffic situation. In a limited traffic- 
related infrastructure scenario, the information from the preceding trip can aid the estimation. Therefore, an 
analysis of the correlation between the preceding trips to the current trip within a 30-minute time frame is 
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conducted initially using (1). The level of variation in one feature due to variation in another is quantified 
through correlation analysis [25] method. The range of the correlation coefficient is -1 to +1, indicating perfect 
negative or positive correlation respectively. 


DI(Cti-vt) *(te-etk)] 
x= 1 
JSC utp*d Ce Hety) (1) 


In (1), t; is the current trip travel time, t, is the preceding trip travel time, pt; and yt, are the means 
of the current and preceding trip travel times and x; is the correlation coefficient. Similarly, the correlation 
coefficient x, of current trip travel speed to that of the preceding trip travel speed is estimated. The coefficients 
estimated are used in developing the proposed model. The Latest Travel Time LTT; ;is estimated using the 
model in (2). 


Qj * TT, j if RS, j < as 


LTT; ; = { 
(245 * TT + 22,5 * TRI || =~ otherwise 


(2) 


The latest travel time LTT; ;is estimated based on the traffic scenario during the previous trip in two traffic 
patterns. 

Case 1: Congestion: If, the running speed of the segment during the previous trip falls below THRS; (the 
threshold speed [26] of the segment j), congestion in the link is presumed. The latest travel time in such a 
situation is calculated by adjusting the preceding trip's travel time TT; ; by the defined correction rate a;. 
Case 2: Congestion-free: If the running speed RS; ; of the segment in the preceding trip is above the threshold 
range, then it is presumed that the link is congestion-free. In such cases, the LTT; ; is estimated based on the 
weighted average of preceding trip travel times TT; ; and TRT;, ; (preceding trip travel time based on speed 
RS,,;). The estimations for @;, 21,;, Z2,;, TRT,,; and THRS; given in (3) to (7). 

The correction rate a; is estimated for the trips of case / type in the historical data. Its value is in the 
range of 0.8 < @ < 2. With the base value of 0.8, it is incremented in the steps of 0.05, up to 2. At each value 
of a, error, is calculated using (3) where the travel time of the current trip is CTT; ; and the preceding trip is 
TT,,; in the archived data. The a@ that yields the minimum error; ; is the a; explored for segment j, and the 
respective a; will be used in (2). 


error; = Yiea|TT,,j *@ — CTT;,;| 08<a<2 (3) 


The threshold for running speed THRS; for each segment j is estimated based on historical data 
through a statistical approach 1.5 * /QR rule as given in (4). Q1 (Quartile 1), is the mean speed value below 
which 25% of data are present, and Q3 (Quartile 3), is the mean speed value below which 75% of them are 
present when they are arranged in ascending order. /QR (InterQuartile Range) is the difference between Q1 
and Q3. Conventionally the constant 1.5 is multiplied by the /QR to estimate the lower and upper bound, the 
data points beyond this are considered outliers. In the current study, only the lower bound is used to predict 
congestion. If the speed of the preceding trip is an outlier (to the lower side), congestion is presumed, and the 
current travel time is increased a; times of the TT, ; 


THRS; = Q1(RS;) — 1.5 * IQR(RS;) "2 


The 2, ; Zz; , are the weights of the travel time and travel speed of the preceding trip respectively. In 
(5) and (6) the x;,,x, are the coefficients estimated using (1). The TRT,; is the travel time for segment j 
estimated based on the running speed RS, ; of the preceding trip using (7), where d; is the distance and RS, ; 
is the preceding trip running speed of the respective segments. The Delay; is the average delay at the 
intersection for segment / and it is added to the travel time estimated through running speed. Delay; is 
determined based on field study and historic location data. 


45 = X¢/(X_ + Xs) (5) 
225 = X5/ (Xp + Xs) (6) 
TRT,,j = seconds(d;/RS,.;) + Delay; (7) 
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3.3. A hybrid model: hybrid travel time estimation model 

The schematic of the implementation process of the hybrid model proposed is presented in Figure 2. 
The proposed HTTEM combines the predictions of the GBRT model and the LTTE model. The HTTEM is 
given in (8), where ETT; ; is the estimated travel time of a trip for segment 7 In this a weighted average of 
predicted travel time PTT; ; by GBRT and the latest travel time LTT; ; by LTTE model, based on the weights 
(W1,;,W2,;). The weights are estimated by applying (10). 


ETT; ; = Wij * PTT;,; + W2,j * LTT; ; (8) 
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Figure 2. The schematic of the proposed model implementation process 


The R-squared (R2) or the coefficient of determination is a statistical parameter that examines the linear 
relationship strength among the variables considered. The general equation for estimating the coefficient of 
determination measure is given in (9). Where y, ; is the actual travel time of trip / for segment / which is the 
dependent variable, x;; is the set of independent variables namely the start time of the trip, day of the week, 
and the LUP, and f(x;,,;) is the function that returns the predicted travel time. wy; is the average of the actual 
travel time in the test data. The R2 value of the ML model is the w, ; in (9), and the formula to estimate wp ; 
is given in (10). 


Z 
Res1- Thy (vi — f(x:;)) [Thais - uy) “ 


j= R? wy = (1-w1,;) (10) 


The step-by-step procedure to dynamically estimate the travel time at the segment level is given in 
Algorithm 1. The steps to estimate the hybrid travel time at the route level are given in Algorithm 2. The 
dynamic travel time at the route level is estimated by summation of the ETT, ; of each segment estimated 
dynamically. 


Algorithm 1: Hybrid travel time estimation procedure at segment level 
Inputs e Start time Start_time,, for segment j 
e Weights w,; and w2; 
e The threshold running speed THRS; for the segment j 
e The correction ratio aj, 4%; and z,; of each segment j 
e Preceding Trip information details: 
a) start time Start_time, ; 
b) Running speed RSx,; 
c) travel time TT,,; 
Output Estimated travel time ETT,; of a trip i for the segment j at Start_time; 
Step 1 At Start_time; predict the travel time PTT,; using the machine model 
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Step 2 Check the traffic pattern if congestion is sensed go to step 4 else estimate the 


Step 3 Estimate the latest travel time and go to step 5 


LTT, ; = 24,5 * TTy,j + 22,5 * TRTy,; 


Step 4 Estimate the latest travel time 


LTT, ; = aj * TT; 


Step 5 Estimate travel time for segment j for the trip i 


ETT, ; = Wy; * PTT;,; + wo; * LTT; ; 


Algorithm 2: Dynamic travel time estimation procedure at route level 
Input Start time Start_time; of the trip i 


Number of segments n 


travel time TRT,; using the RS&,,;, the running speed of the preceding trip 


Preceding Trip information details of n segments along the route: 
a) start time Start_time,, j, b) travel time TT,,;, c) Running speed RS,,; 
Output Estimated travel time for the route RTT; with start time Start_time; 


Step 1 Initialize the start time of the initial segment x= Start_time;, 


uy 


where j=1 


Step 2. Fore cast PTT,; for segment j using the GBRT model with the start time as x 


Step 3 Estimate the latest travel time 


Ty es ee «TT, j +22; *TRT pj 


aj * TT, if RS, j < ae 


otherwise 


Step 4 Using the proposed Hybrid model calculate the travel time of the segment j 


ETT, 5 = Wy; * PTT; + Wo; * LTT, j 


ETT, = ETT; 


Step 5 Estimate the start time of the next segment 


x =X + ETT,; 


j=jtl 


repeat steps 2 to 5 for each segment 
Step 6 Estimate the route travel time at the start time 


k 


RTT sti = ETT, 
jz 


4. RESULTS AND DISCUSSIONS 


For estimating the travel time using the proposed HTTEM, the correlation coefficients, the threshold 
running speed, the correction rate, and the R2 of the GBRT predictions are estimated using the models 
presented in the previous section, and the numerical results are presented in Table 3. A positive correlation is 
observed between the preceding to current trip travel times. The threshold speeds for segments vary from 
12-25 km/hr indicating the impact of the LUP on travel time. The correction rate varies between 0.95 to 1.1. 


Table 3: Coefficients and constants estimated in the study 


Xt Xs (21 j, 2Z2;) (Wj, Wa;) aj THRS; 
Segment_l 0.52 0.65 (0.46,0.54) (0.69,0.31) 0.95. 12 
Segment_2 0.69 0.52  (0.57,0.43) (0.63,0.33) 1.05. 15 
Segment_3 0.50 0.46 (0.52,0.48) (0.65,0.35) 11 18 
Segment_4 0.49 0.58 _(0.44,0.56) _(0.62,0.38) 1.0 25 


4.1. Predictions at the route and segment levels 

The predictions of the GBRT, LTTEM, and the proposed HTTEM against the actual travel time are 
presented in this section. Algorithm | is implemented and the travel times are estimated for each segment. The 
model predictions for a few samples at various start times of the day of all segments are presented in Figure 3. 
It is observed that the travel time estimates of the HTTEM are close to the actual as compared to LTTEM and 
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GBRT models emphasizing the suitability of a hybrid model. Algorithm 2 is implemented to estimate the travel 
time at the route level. The predictions of each model at the route level are presented in Figure 4. It is observed 
from the plots that, the predictions of the proposed HTTEM are superior as compared to GBRT and LTTEM 
highlighting the adaptability of the hybrid model at various times of the day. 
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Figure 3. The actual and predicted travel times by GBRT and the proposed models for all segments 
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Figure 4. The actual and predicted travel times by GBRT and the proposed HTTEM at the route level 


4.2. Model performance 

The error in the predictions of the GBRT model, LTTEM, and HTTEM are estimated based on the 
most commonly used regression error metrics [27] such as median absolute deviation (MAD), mean squared 
error (MSE), root mean squared error (RMSE), mean absolute error (MAE) and mean absolute percentage error 
(MAPE). The errors are summarized in Table 4. It is observed from the errors that the HTTEM outperforms 
the LTTEM and GBRT models. It is also observed that percentage errors at coarse spatial aggregation i.e., 
route level is better than at segment level which is a fine aggregate. This indicates that at a higher spatial 
aggregation level the variations are normalized, and the predictions are better. There is also a possibility that 
the variations and the delays experienced at one segment by the buses are compensated by the driver during 
traversing [11] the other segments. For applications like route travel time predictions and timetable [28] 
generation, predicting travel times at route level spatial aggregation is needed and the proposed model is 
recommended. 

At fine spatial aggregation levels, the impact of the variations is apparent, but this gives deeper 
insights into the travel times. For applications like bus arrival time prediction at bus stops and other passenger 
information systems, the predictions at the lower spatial levels are preferred for making precise predictions. 
Further, to optimize the travel times of a route, if the treatments or solutions are provided to individual segments 
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will have a better impact as the spatial properties are different for each segment. The predictions made at the 
segment level can be extended to bus stop levels in the future. Also, as a supplement to the error metrics to 
demonstrate the performance of the models, two new error scales were used. 

— Percentage of trips with less than 60 seconds of error. 

— Percentage of trips with less than 30 seconds of error. 

The percentage of trips in both error scales is presented in the later columns of Table 4. The results 
depict that when the HTTEM is used on average 92.5% of trips have errors of less than a minute (60 seconds) 
at segment level and 66% at route level which are better as compared to GBRT with 84% at segments and 60% 
at the route level and LTTEM with 86% at the segment and 64% at the route level. Similarly, the errors of the 
HTTEM are less than 30 seconds during 78.5% of the trips at the segment level and 39% at the route level, 
respectively, whereas it is 61.5% and 30% for the GBRT model and 68% and 32% for LTTEM. From all the 
comparisons made in this section, it is evident that the proposed HTTEM outperforms the LTTEM and GBRT 
models. This also emphasizes the importance of the input data used for modeling. The proposed HTTEM and 
LTTEM perform better as compared to the GBRT model, as these models are developed based on recent trip 
information and historic data whereas the GBRT model uses only historic data. This highlights the need for 
dynamic data for short-term travel time prediction. 


Table 4: Prediction errors comparison of the three models 
Spatial entity Error metrics Error in seconds and percentage = Error scales _ Percentage of trips in error scales 


GBRT _LTTEM _HTTEM GBRT _LTTEM _ HTTEM 

Route MAD 49.0 47.0 48.0 < 60 seconds 60% 64% 66% 
MAE 64.0 52.0 52.0 

RMSE 84.0 72.0 60.0 < 30 seconds 30% 32% 39% 
MAPE 5.1% 4.8% 4.0% 

Segment-1 MAD 16.5 17.7 13.7 < 60 seconds 82% 81% 88% 
MAE 22.9 18.8 18.4 

RMSE 31.6 23.7 24.3 < 30 seconds 61% 60% 66% 
MAPE 5.4% 4.6% 4.5% 

Segment-2 MAD 27.4 24.1 19.5 < 60 seconds 91% 89% 93% 
MAE 32.4 24.5 21.6 

RMSE 39.4 28.0 26.3 < 30 seconds 60% 58% 63% 
MAPE 13.5% 10.1% 9.1% 

Segment-3 MAD 13.4 17.0 13.7 < 60 seconds 73% 80% 93% 
MAE 14.9 20.2 13.7 

RMSE 18.4 25.8 16.0 < 30 seconds 62% 85% 91% 
MAPE 5.2% 7.1% 4.9% 

Segment-4 MAD 16.5 19.0 14.6 < 60 seconds 89% 92% 96% 
MAE 18.4 20.1 15.6 

RMSE 21.8 23.5 19.0 < 30 seconds 63% 10% 73% 
MAPE 8.0% 8.7% 6.9% 


The authors call attention to the need for real-time road traffic information through sensors, 
crowdsourced data of commuters [29], passenger counters, and bus passenger demand information to further 
optimize the dynamic travel time predictions. Multiple data sources can improve the model performance but 
integrating multiple data sources for real-time applications is again a challenge. Setting the infrastructure for 
dynamic traffic flow information across the entire road network is also economically not viable in a few 
countries. The crowdsourced data is a promising and economically viable supplementary data set. The current 
work is one such case where the regular bus trips are used as the probe vehicle trips and the details of the 
preceding trip as dynamic traffic information on the roads. Hence for such urban areas and tier-2 cities, the 
proposed model can be a feasible option with considerable performance to predict travel time. 


5. CONCLUSION 

In this work, a HTTEM for dynamic travel time prediction is proposed. Sample trips of Tumakuru 
City Service, India during weekdays between mornings 7:00 to evening 20:00 are used for modeling. The 
location data are processed to trip aggregates at two spatial levels; route and segment, based on land use 
patterns. The GBRT is trained with historic data, to predict the travel time for the current start time of a trip at 
both route and segment levels. A LTTEM is developed based on preceding trip details and historical data. The 
proposed HTTEM dynamically adjusts the results of the GBRT model and the LTTEM. The predictions made 
by the model were assessed using error metrics MAD, MSE, RMSE, and MAPE and two additional error scales 
illustrating the percentage of trips that estimate the short-time travel time with less than 60 and 30 seconds of 
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error. The performance of the proposed HTTEM demonstrated better accuracy at the route and segment level 
as compared to the GBRT and LTTEM. The authors also emphasize supplementing the location data with 
crowdsourced data sources as an economically viable option to optimize the dynamic travel time predictions 
in the future. Overall, it is concluded that, in urban areas and tier-2 cities where there are limited traffic-related 
data sources, the proposed model that uses the preceding buses as probe vehicles can predict decent results. 
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